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Abstract. We introduce a general algebraic setting for describing linear 

boundary problems in a symbolic computation context, with emphasis 

Jh ^ on the case of partial differential equations. The general setting is then 

applied to the Cauchy problem for completely reducible partial differ- 
ential equations with constant coefficients. While we concentrate on the 
1^^, , theoretical features in this paper, the underlying operator ring is imple- 

rvj ' mented and provides a sufficient basis for all methods presented here. 



^ ■ 1 Introduction 

j/2 ■ A symbolic framework for boundary problems was built up in |11I13| for linear 

O I ordinary differential equations (LODEs); see also |15I6I7| for more recent devel- 

opments. One of our long-term goals is to extend this to boundary problems for 
linear partial differential equations (LPDEs). Since this is a daunting task in 
^ . full generality, we want to tackle it in stages of increasing generality. In the first 

^^ ' instance, we restrict ourselves to constant coefficients, where the theory is quite 

well-developed [3 . Within this class we distinguish the following three stages: 

C ] ' 1. The simplest is the Cauchy problem for completely reducible operators. 

^^ , 2. The next stage will be the Cauchy problem for general hyperbolic LPDEs. 

^D ' 3. After that we plan to study boundary problems for elliptic/parabolic LPDEs. 

m ■ 

In this paper we treat the first case (Section 2]). But before that we build 

up a general algebraic framework (Sections [5] and [3]) that allows a symbolic 

description for all boundary problems (LPDEs/LODEs, scalar/system, homoge- 

/\ . neous/inhomogeneous, elliptic/hyperbolic/parabolic). Using these concepts and 

tools we develop a general solution strategy for the Cauchy problem in the 
case (1). See the Conclusion for some thoughts about the next two steps. 
The passage from LODEs to LPDEs was addressed at two earlier occasions: 

— An abstract theory of boundary problems was developed in [5], including 
LODEs and LPDEs as well as linear systems of these. The concepts and re- 
sults of Sections [2] and |3] are built on this foundation, adding crucial concepts 
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whose full scope is only appreciated in the LPDE setting: boundary data, 
semi-homogeneous problem, state operator. 
— An algebraic language for multivariate differential and integral operators was 
introduced in Section 4 of 14 , with a prototype implementation described 
in Section 5 of the same paper. This language is generalized in the FIDOS 
algebra of Section |4l and it is also implemented in a Mathematica package. 

In this paper we will not describe the current state of the implementation 
(mainly because of space limitations). Let us thus say a few words about this 
here. A complete reimplementation of the FIDOS package described in [T3] is 
under way. The new package is called OFIDO (Ordinary and Fartial Integro- 
Differential Operators), and it is implemented as a standalone Mathematica 
package unlike its predecessor, which was incorporated into the THEOREMA 
system. In fact, our reimplementation reflects several important design princi- 
ples of THEOREMA, emphasizing the use of functors and a strong support for 
modern two-dimensional (user-controllable) parsing rules. We have called this 
programming paradigm FUNFRO, first presented at the Mathematica sympo- 
sium [12]. The last current stable version of the (prototype) package can be 
found at |http : //www . kent . ac . uk/smsas/personal/mgr/index . htmH 

At the time of writing, the ring of ordinary integro- differential operators is 
completed and the ring of partial integro- differential operators is close to com- 
pletion (for two independent variables). Compared to [13], the new FIDOS ring 
contains several crucial new rewrite rules (instances of the substitution rule for 
resolving multiple integrals). Our conjecture is that the new rewrite system is 
noetherian and confluent but this issue will be analyzed at another occasion. 

Notation. The algebra of m x n matrices over a field K is written as K™, 
where m = 1 or n = 1 is omitted. Thus we identify K"^ = K ® ■ ■ ■ ® K with the 
space of column vectors and Kn = {K"^)* with the space of row vectors. More 
generally, we have K™ = K" -^ K"\ 

2 An Algebraic Language for Boundary Data 

As mentioned in the Introduction, we follow the abstract setting developed in [9] . 
We will motivate and recapitulate some key concepts here, but for a fuller treat- 
ment of these issues we must refer the reader to 19, and its references. 

Let us recall the notion of boundary problem. Fix vector spaces J- and Q over 
a common ground field K of characteristic zero (for avoiding trivialities one may 
assume J-' and Q to be infinite-dimensional). Then a boundary problem (T,B) 
consists of an epimorphism T: T —^ Q and a subspace B C T* that is orthogo- 
nally closed in the sense defined below. We call T the differential operator and B 
the boundary space. 

Similar to the correspondence of ideals/varieties in algebraic geometry, we 
make use of the following Galois connection [9l A. 11]. If ^ is any subspace of the 
space F, its orthogonal A^ < T* is defined as {^p ^ T* \ (p{a) = for all a G A}. 
Dually, for a subspace B of the dual space J-"* , the orthogonal B-^ < J-" is defined 
hy {f G T\ I3{f) = for all l3 e Bj.liwe think of J" as "functions" and of T* as 



"boundary conditions" , then A^ is the space of vahd conditions (the boundary 
conditions satisfied by the given functions) while Z?^ is the space of admissible 
functions (the functions satisfying the given conditions). 

Naturally, a subspace of either 5 of J^ or J-'* is called orthogonally closed if 
<S^^ = S. But while any subspace of J-" itself is always orthogonally closed, this 
is far from being the case of the subspaces of the dual T* . Hence the condition on 
boundary spaces B to be orthogonally closed is in general not trivial. However, 
if B is finite-dimensional as in boundary problems for LODEs (as in Example [1] 
below), then it is automatically orthogonally closed. For LPDEs, the condition 
of orthogonal closure is important; see Example [5] for an intuitive explanation. 

In [11113) and also in the abstract setting of |1] we have only considered 
what is sometimes called the semi-inhomgeneous boundary problem [16|, more 
precisely the semi-inhomogeneous incarnation of {T,B); see Definition [7] for the 
full picture. This means we are given a forcing function f ^ Q and we search for 
a solution u & J- with 



P{u) 



/, 



(/3 e B). 



(1) 



In other words, u satisfies the inhomogeneous "differential equation" Tu — f 
and the homogeneous "boundary conditions" /3(m) = given in B. 

A boundary problem which admits a unique solution u & T for every forcing 
function / G ^ is called regular. In terms of the spaces, this condition can be 
expressed equivalently by requiring that Ker T + B'^ = J-\ see [2| for further 
details. In this paper we shall deal exclusively with regular boundary problems. 
For singular boundary problems we refer the reader to [5] and [S]. 

For a regular boundary problem, one has a linear operator G : ^ — > J^ send- 
ing / to u is known as the Green's operator of the boundary problem (T,B). 
From the above we see that G is characterized by TG — 1 and Im G = B^ . 

Example 1. A classical example of this notion is the two-point boundary prob- 
lem. As a typical case, consider the simplified model of stationary heat conduction 
described by 



u" = /, 

u(0) = m(1) = 0. 



Here we can choose T ^ Q = C°° (R) for the function space such that the differ- 
ential operator is given hy T — D"^ : C°° (R) -^ C°° (R) and the boundary space 
by the two-dimensional subspace of C°°(R)* spanned by the linear function- 
als L: u I— !> u(0) and R: u^^ u{l) for evaluation on the left and right endpoint. 
In the sequel we shall write B = [L,R], employing an important generalization 
for LPDEs, described in the next camples. We can express its Green's operator 
in the language of integro-differential operators as explained in |13j . 

Example 2. As a typical counterpart in the world of LPDEs, consider the equa- 
tion for waves in an inhomogeneous medium, described by 



Uxx -Utt = f{x,t) 

u{x, 0) = ut(x, 0) = u(0, t) ^u{l,t) = 



in one space dimension. In this case we choose J- = Q — C"^(R x IR+); again one 
could choose much larger spaces of functions (or distributions) in analysis and in 
the applications. Here the differential operator is T = Dxx — Du '■ C"^(]RxR"'") — )■ 
C"^(]R X R+) while the boundary space B is the orthogonal closure of the linear 
span of the families of functional Px,lx {x e IR) and Kt,Xt;{t G IR+) defined 
by I3x{u) = u{x,0),jx{u) = ut{x,0) and Kf(u) = u{0,t),Xt{u) = u{l,t). Using 
the notation [. . . ] for denoting the orthogonal closure of the linear span, we can 
thus write B = [/3a;,7x, k^, At | x £ M.,t G R+] for the boundary space under 
consideration. 

The point of the orthogonal closure is that the given conditions imply other 
conditions not in their span, for example Ux{l/2,0) = or j _^u{0,t) dr = 0. 
Rather than being linear consequences, these two examples are differential and 
integral consequences. (Of course the full boundary space also contains many 
functionals without a natural analytic interpretation.) 

In the problems above, the differential equation is inhomogeneous while the 
boundary conditions are homogeneous. A semi-homogeneous boundary problem 
is the opposite, combining a homogeneous differential equation with inhomoge- 
neous boundary conditions. While this is a simple task for LODEs (as always we 
assume that the fundamental system is available to us in some form!), it is usu- 
ally a nontrivial problem for LPDEs (even when they have constant coefhcients). 
We will give the formal definition of a semi-inhomogeneous boundary problem 
in the next section (Definition [7|) . Here it suffices to consider an example for 
developing the necessary auxiliary notions. 

Example 3. The Cauchy problem for the wave equation in one dimension is 



Uxx - Utt = 0, 

u{x, 0) = f{x), ut{x, 0) = g{x). 



Being a hyperbolic problem, we could use rather general function spaces for 
the "boundary data" /, g. For reasons of uniformity we will nevertheless restrict 
ourselves here to the analytic setting, so assume f,g € C^(R). Note that the 
association of u to (/, g) is again a linear operator mapping (two univariate) 
functions to a (bivariate) function; we will come back to this point in Definition[7l 

Going back to the abstract setting, one is tempted to define the notion 
of boundary data as some kind of functions depending on "fewer" variables. 
But the problem with this approach is that — abstractly speaking — we are not 
dealing with any functions depending on any number of variables (but see be- 
low). Moreover, the inhomogeneous boundary conditions in the form u{0,x) = 
f{x),ut{0,x) — g{x) are basis-dependent while the whole point of the abstract 
theory is to provide a basis-free description (which leads to an elegant setting 
for describing composition and factorization of abstract boundary problems); see 
Proposition [HI We shall therefore develop a basis-independent notion of bound- 
ary data (we can go back to the traditional description by choosing a basis). 



We define first tlie trace map tre: J-" — >■ B* as sending / G J-" to the func- 
tional /3 i-^ /?(/)• In Example 12] this would map the function u{x, t) to its position 
and velocity values on R x {0}. We call trc(/) the trace of / and write it as /*. 
Moreover, we denote the image of the map trc by B' and refer to its elements as 
boundary data. Note that B' is usually much smaller than the full dual B* since a 
continuous function (let alone an analytic one) cannot assume arbitrary values. 
(This situation is vaguely reminiscent of the algebraic and continuous dual of a 
topological vector space.) 

Since by definition the trace map is surjective from 7^ to B', it has some 
right inverse B^ : B' — > J-". We refer to B^ as an interpolator for B since it 
constructs a "function" / — B^{B) e T from given boundary values B E B' 
such that /?(/) = B{P). Of course, the choice of / is usually far from being 
unique. Apart from its use for describing boundary data (see at the end of this 
section), the notion of interpolator will turn out to be useful for solving the 
semi-homogeneous boundary problem (see Proposition [TU]). 

Let us now describe how to relate these abstract notions to the usual setting 
of initial and boundary values problems as they actually in analysis: essentially 
by choosing a basis. However, we have to be a bit careful since we must deal 
with the orthogonal closure. 

Definition 4. li B < T* is any orthogonally closed subspace, we call a fam- 
ily (/3i I « G /) a boundary basis ii B = [I3i \ i £ I], meaning B is the orthogonal 
closure of the span of the /3i . 

Note that a boundary basis is typically smaller than a X-linear basis of B. 
All traditional boundary problems are given in terms of such a boundary basis. 
In Example El the boundary basis could be spelled out by using / = ]R l+) R 
with /3(a;^o)(w) — u,{x,0) and /3(2,_i)(m) — Ut{x,0). Relative to a boundary ba- 
sis {Pi \ i G I), we call / — /3i{f)i(zi G K^ the boundary values oi f € F. As we 
can see from the next proposition, we may think of the trace as a basis-free de- 
scription of boundary values. Conversely, one can always extract from any given 
boundary data B € B' the boundary values B{/3i)i^j as its coordinates relative 
to the boundary basis (/3i). 

Lemma 5. Let B < J-* be a boundary space with boundary basis {Pi |g /). // 
for any B,B G B' one has B{Pi)i^i — B{Pi)i^i then also B — B. In particular, 
for any f G J-, the trace f* depends only on the boundary values f{Pi)iizj. 

Proof. Since B' is the image under the trace map, we have B — f* and B = f* 
for some f,f£T. So assume f{Pi) — f{Pi) for alH G / and thus /?(/ — /) = 
for all /3 G B by the definition of orthogonal closure. Then we have trc(/ — /) == 
and thus /* = /*. D 

The analytic interpretation of this proposition is clear in concrete cases like 
Example[2J Once the values u{x,Q),ut{x,Q) and u(0,i), u(l, t) are fixed, all dif- 
ferential and integral consequences, as in the above examples Ma;(l/2,0) = 
or Jj^^^u(0,t) (ir, are likewise fixed. It is therefore natural that an interpolator 



need only consider the boundary values rather than the full trace information. 
This is the contents of the next lemma. 

Lemma 6. Let B < J-* he a boundary space with boundary basis {j3i |g /) 
and write fi = f{(3i)i^i e K^ for the boundary values of any / G J-" and Bj 
for the K-subspace of K^ generated by all boundary values fi. Then any linear 
map J: Bj —^ J^ with J{fi)i — fi induces a unique interpolator B^ : B' ^ F 
defined by B i-^ J{B(j3i)i^i). 

Proof. We must show that B^ : Z?' — s> 7^ is a right inverse of trc: T ^r B' . So 
for arbitrary B =G B' we must show J{B{(3i)i(zi)* — B. By the definition of B' 
we can write B = f* for some f G J^. Since f*{(3i) = (3i{f), we are left to 
prove J{fi)* — /*. Using Lemma[Sl it suffices to prove that J{fi)i — fi, which 
is true by hypothesis. D 

As noted above, we can always extract the boundary values B{[3i)i^j £ K^ of 
some boundary data B £ B' relative to fixed basis (/3i) of B. However, since one 
normally has got only the boundary values (coming from some function) , where 
does the corresponding B € B' come from? By definition, it has to assign values 
to all ^ G B, not only to the (3i making up the boundary basis. As suggested 
by the above lemmata, for actual computations those additional values will be 
irrelevant. Nevertheless, it gives a feeling of confidence to provide these values: 
If B^ is any interpolator, we have B{(3) — f3{B'^{Bi)i^i). This follows immedi- 
ately from the fact that B"^ is a right inverse of the trace map and that it depends 
only on the boundary values {Bi)i^i by Lemma [B] In the analysis setting this 
means we interpolate the given boundary value and then do with the resulting 
function whatever is desired (like derivatives and integrals in Example ^ . 

3 Green's Operators for Signals and States 

Using the notion of boundary data developed in the previous section, we can 
now give the formal definition of the semi-homogeneous boundary problem. In 
fact, we can distinguish three different incarnations of a "boundary problem" 
(as we assume regularity, the fully homogeneous problem is of course trivial) . 

Definition 7. Let {T,B) be a regular boundary problem with T: T ^ Q and 
boundary space B C T* . Then we distinguish the following problems: 

Given (/, B)£g®B', Given /eg. Given B £ B' , 

find u £ F with find u £ T with find u £ T with 



Tu = f, 

l3{u) = B{l3){l3£B). 



Tu = f, 

P{u) = (/3 G B). 



Tu^O, 

/3(u) = B{I3) (13 £ B). 



They are, respectively, called the fully inhomogeneous, the semi-inhomogeneous 
and the semi-homogeneous boundary problem for (T, B). The corresponding lin- 
ear operators will be written a.s F : Q (B B' —>■ J^ , {f,B) i— )• u and G: Q ^ F, 
f 1-^ u and H: B' ^ F, B h^ u. 



Lemma 8. Each of the three problems in Definition[^has a unique solution for 
the respective input data, so the operators F, G, H are well-defined. 

Proof. Employing the usual superposition principle F — G ® H , we can restrict 
ourselves to the semi-inhomogeneous and the semi-homogeneous boundary prob- 
lem. For the former, the existence and uniqueness is verified in W, §2]. For the 
semi-homogeneous problem, existence is seen as follows: Since the given bound- 
ary data B E B' can be written as B — f* by the definition of B' , we have to 
find u € J- such that Tu = and /3(u) = j3{f) for all fS € B. Using the decompo- 
sition / = fo+fi corresponding to the direct sum Ker T+B-^ = T, we set u = /q. 
Then Tu = is clear, and we check further that /3(u) = /?(/)- /?(/i) = /3(/). For 
uniqueness, it suffices to prove that Tu — 0,u* = has only the trivial solution, 
which is clear since the sum KerT -j- B-^ = J^ is direct. D 

The terminology for the operators F, G, H is not uniform in the literature. In 
the past, we have only considered G and called it the "Green's operator" acting 
on a "forcing function" /. While this is in good keeping with the engineering 
tradition and large parts of the standard mathematical culture |16j , it is difficult 
to combine with suitable terminology for F and H . In this paper, we shall follow 
the systems theory jargon [5 ^^'^ refer to F as the (full) transfer operator, 
to G as the (zero-state) signal transfer operator or briefly signal operator, and 
to H as the (zero-signal) state transfer operator or briefly state operator. This 
terminology reflects the common view of forcing functions / G J-" as "signals" 
and boundary data B Q B' a,s (initial) "states" . 

One of the advantages of the abstract formulation is that it allows us to 
describe the product of boundary problems in a succinct, basis-free manner (and 
it includes LODEs and LPDEs as well as systems of these). The composite 
boundary problem can then be solved, both in its semi-inhomogneous and its 
semi-homogeneous incarnation (the latter is presented here for the first time). 

Proposition 9. Define the product of two boundary problems (T, B) and (T, B) 
with F^g ^n andBCg*, BCF by 

{T,B){f,B) = {Tf,Bf + B). 

Then (T,B){T,B) is regular if both factors are. In that case, if (T,B), {T,B) 
have, respectively, the signal operators G, G and the state operators H, H, 
then {T, B){T, B) has the signal operator GG and the state operator (BT + B)' — >■ 
J" acting by B + B ^ GH{Bf*) + H{B). 

Proof. The preservation of regularity and the relation for the signal operators 
is proved in Propositions. 2 of [S]- For the statement about the composite state 
operator, note first that the sum BT+B is direct by [H (3.2)]; hence the definition 
is consistent. Now let B -(- B G {BT + B)' be arbitrary boundary data and 
set u = GH{Bf*) + H{B). Then TTu = TH{Bf*) = since TG = 1 and 
both TH and TH vanish by the definition of state operator. Hence the differential 
equation of the composite semi-homogeneous boundary problem is satisfied. 



It remains to check the boundary conditions /3f{u) = B{(3T) for /3 G i3 
and /3(u) = i?(/3) for JS e B. For the first, we use TG = 1 and TH^O again 
to compute /3f(u) = l3H{Bf*). Since f *B e i3' and H is the state operator 
for iT,B), we obtain l3H{Bf*) = {Bf*){(i) = ^B{l3f) as required. For the 
second set of boundary conditions, we use that /3G = since G is the signal 
operator with homogeneous boundary conditions P G B. Hence (3{u) — (3H{B), 
and now the claim follows because H is the state operator for the inhomogeneous 
boundary conditions /S € B. D 

As detailed in jlll9| , the computation of the signal operator G can be decom- 
posed in two parts: (1) Finding a right inverse T* of the differential operator T, 
which involves only the differential equation without boundary conditions (so 
we may replace the boundary by intial conditions, thus having again a unique 
solution: this is the so-called fundamental right inverse). (2) Determining the 
projector onto the homogeneous solution space along the space of functions ad- 
missible for the given boundary conditions — the projector "twists" the solutions 
coming from the right inverse into satisfying the boundary conditions. An anal- 
ogous result holds for the computation of the state operator H if we replace the 
right inverse T* of T by the interpolator B^ for the boundary space B. 

Proposition 10. Let (T, B) be regular with operators F, G, H as in Definition^ 
Then we have G = (1 - F) T* and H = PB^ , hence F = (1 - P) T* © PB^ for 
the transfer operator. Here T* : Q ^ J- is any right inverse of the differential 
operator T : J^ ^f Q and B^ : B' ^ T any interpolator for B while P : J- ^ J- is 
the projector determined by Im P = Ker T and Ker P = B^ . 

Proof The formula for G is given in [2 (2.3)]. For proving H = PB^ , let B ^ B' 
be arbitrary and set u — PB^{B). Then Tu = follows since ImP = KcrT. 
Furthermore, for every jS £ B we have 

/3(w) = ^(B*P) - /3((1 - P)B'^B) = B{P) - 

by the definition of ;B* and Im(l-P) = KcrP = B^ . This means that u = H{B) 
satisfies the boundary conditions, so H solves the semi-homogeneous boundary 
problem for (T,i3). D 

If T is a completely reducible differential operator with constant coeffi- 
cients in (D , the determination of T^ reduces to solving an inhomogeneous first- 
order equation with constant coefficients — which is of course straightforward 
(Lemma [TU). Also the determination of the interpolator B^ turns out to be easy 
for a Cauchy problem since it is essentially given by the corresponding Taylor 
polynomial ^. Hence it remains to find some means for computing the kernel 
projector P for a boundary problem (T, B) . 

In the case of a LODE of order n, the method for computing P given in 
the proof of Theorem 26 of [13] and in Section 6 of [5] is essentially a Gaussian 
elimination on the so-called evaluation matrix I3{u) = [I3i{uj)]ij G X"^" formed 
by evaluating the i-th boundary condition l3i on the j-th fundamental solution Uj . 



So here we assume ui, . . . , u„ is a basis of KerT and /3i, . . . , /3„ a basis of B. 
Unfortunately, this is not a very intuitive description of P, and it is not evident 
how to generalize it to the LPDE case. We have to gain a more conceptual 
perspective at /3(m) for making the generalization transparent. 

Let us write Ev : B Ker T -^ K for the bilinear operation of evaluation 
{13, u) H' /3(m). Choosing bases /3i,...,/3„ for B and ui,...,'u„ for KerT, the 
coordinate matrix of Ev is clearly /3(m). By the usual technique of dualization, 
we can also think of Ev as the map B : Ker T ^- B* that sends u e Ker T to 
the functional /? i-^ P{u). But this map is nothing else than the restriction of 
the trace map trc : T -^ B' to Ker T C J-". It is easy to check that the restricted 
trace is bijective and that its inverse gives rise to the projector. 

Proposition 11. Let (T, B) be a regular boundary problem with E: KerT — > B' 
being the restricted trace map. Then E is bijective with the state operator H as 
its inverse, and P — Ho trc is the projector with ImP = KerT and KerP = B'^ . 

Proof. Given any boundary data B G B' , we know that u = H{B) satis- 
fies the inhomogeneous boundary conditions (3{u) = B{f3) so that u* = B. 
Hence ^: KerT ^ B' is a left inverse oi H : B' -> KerT C J". But it is 
also a right inverse because u £ KerT is a solution of the semi-homogeneous 
boundary problem with boundary data B = u* ^ B' and thus must coincide 
with H{B) E KerT by the uniqueness of solutions. 

For the projector, observe first that trc: T —^ B' is a left inverse oi H: B' ^ 
T . As is well known 9, (A. 16)], this implies at once that P — ff otrc: T ^ T \s 
a projector with Im P = Ini H and Ker P — Ker(trc) . It remains to show Im 11 = 
KerT and Ker(trc) — B^ . The second identity follows from the definition of trc. 
For the first identity, the inclusion from left to right follows from the definition 
of -ff, the reverse inclusion by writing u ~ H{Eu) for u G KerT. D 

We observe that the formula P{u) ~ H{u*) has a very natural interpretation: 
The kernel projector picks up the boundary data u* of an arbitrary function u G 
J^ and then constructs the required kernel element H{u*) S KerT by solving 
the semi- homogeneous boundary problem with boundary data u* . In the LODE 
case, the relation P — H o trc reduces to the aforementioned formulae (see [5] 
after Proposition6.1) after choosing bases ui, . . . , u„ for KerT and /3i, . . . , /J„ 
for B. Apart from its conceptual clarity, the advantage of PropositiorfTTl is that 
it can also be used in the LDPE case (see after Lemma [Til) . 

4 The Cauchy Problem for Analytic Functions 

At this point we switch from the abstract setting of Sections [5] and [3] to the 
concrete setting of analytic functions. Note that we are dealing with complex- 
valued functions of real arguments. This means the ground field is K — (D, and J- 
is the integro-diflFerential algebra of entire functions restricted to real arguments. 
More precisely, we shall employ the following conventions for easing the bur- 
den of book-keeping: As elements of J- we take all holomorphic functions R" — >■ 



(D for any n e N, including the constant functions u e (D for n = 0. In other 
words, J-" is a direct hmit of algebras. Moreover, we have derivations _D„ and 
integrals A„ for all n > 0, namely D„{u) = du/dxn and 



A„(u)= / u{...,^,...)d^, 
Jo 

where ^ occurs at the n-th position. Clearly, we have then integro-differential 
algebras {J-,D„,An) for every n > 0. In fact, T has the structure of a hierar- 
chical integro-differential algebra. This notion will be made precise at another 
occasion; for the moment it suffices to make the following observations. If N® is 
the sublattice of the powerset ■p(N''') that consists of finite sets a — {ai, . . . , ak} 
and the full set N+ = {1,2,3,...}, we define for a G N® the subalgebras 

J'a = {f eJ'lDif = for all i i a}, 

consisting of the functions depending (at most) on x^j, . . . ,a;Q,^. Then (J^q, C) 
is a sublattice of (7^, C) that is isomorphic to the lattice (N®, C). The bottom 
element is of course J^0 = (D, the top element T^+ = F. We write J^„ as an 
abbreviation for J-{i,,,,^n}- 

As in the earlier paper [T3], we add to this algebraic structure all linear 
substitution operators. In accordance with the above hierarchical structure, we 
use the ring (D* of row and column finite matrices with complex entries|j This 
means any M G (D* can actually be seen as a finite matrix M e (D™ with m 
rows and n columns, extended by zero rows and columns. As usual, we iden- 
tify Af e €™ with the linear map M: C" -> C™, yielding the substitution 
operator M* : T^ — > T^ defined by u{x) i-> u{Mx). 

We write J^[D, A] for the FIDOS algebra generated over (D by the opera- 
tors Dm An {n > 0), the substitutions M* induced by M G (D* and the ex- 
ponential basis polynomials x^e^^ G J^. Here x denotes the arguments x = 
{xi, . . . , Xn) for any n > 0, with exponents a = (ai, . . . , a„) G (D and frequen- 
cies A = (Ai, . . . , A„) G (D. Obviously, J-[D, A] acts on T, with D = (Di, D2, . . . ) 
acting as Di, D2, ■ • ■ and A = {Ai,A2, . . . ) as Ai,A2, . . . , similar to the uni- 
variate case in the older notation of [11]. Here we avoid the notation J for the 
integrals since the powers J might be mistaken as integrals with upper bound n. 

The algebra J^[D, A] can be described by a rewrite system (FIDOS = partial 
integro-differential operator system), analogous to the one given in [T3]. We 
will present this system in more detail — in particular proofs of termination and 
confluence — at another occasion. 

Since in this paper we restrict ourselves to the analytic setting, we can appeal 
to the well-known Cauchy-Kovalevskaya theorem [101 Thm. 2.22] for ensuring 



^ The usage of complex substitutions in functions of a real argument may sound strange 
at first. But an analytic function on R" is of course also analytic on C" with values 
in C so there is no problem with this view. For example, the substitution (l,i)* 
sends f{x) = e^ G J^i to f{x + iy) — e^ cosy 4- ie^ siny G J-2. Moreover, complex 
substitutions are indispensible for specifying the general solution of elliptic equations 
like the Laplace equation. 
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the existence and uniqueness of the solution of the Cauchy problem. While the 
theorem in its usual form yields only local results, there is also a global version [H 
Thm. 7.4] that provides a good foundation for our current purposesjj Since this 
form of the theorem is not widely known, we repeat the statement here. 

As usual, we designate one lead variable t, writing the other ones xi,X2, ■ ■ ■ 
as before. Note that in applications t is not necessarily time. The apparently 
special form of the differential equation Tu = implies no loss of generality: 
Whenever T e <E[Z3] is a differential operator of order m, the change of vari- 
ables t — t,Xi = Xi + t leads to an equation of the required form. 

Theorem 12 (Global Cauchy-Kovalevskaya) . Let T e (D[A,£'i, . . . ,D„] 
he a differential operator in Caucy-Kovalevskaya form with respect to t, mean- 
ing T = _D™ -\-T with deg(T, t) < m and deg(T) < m. Then the Cauchy problem 



Tu = 



(2) 



Dl u{0,xi,...,Xn) = f^ixl,...,Xn) for i:^ l,...,m J 
has a unique solution u G J^n^i for given (/i, . . . , /„) £ J^™. 

In the abstract language of Sections [2] and [3] this is the semi-homogeneous 
boundary problem (T, B) with boundary space 

B = [Lq^^DI I i = 0, . . . , m - 1 and ^ G R™], 

where the evaluation w(t, xi, . . . , Xn) >— >■ w(0, ^i, . . . , ^„) is written as the substi- 
tution Lo,4 = diag(0,^i, . . . ,^m)* denotes . Hence the solution of ^ is given 
by the state operator (/i,...,/,„) G J^™ i— > u if we identify the boundary 
data B £ B' with its coordinate representation (/i, . . . , /„j) G J-"™ relative to the 
above boundary basis {Lq,^DI). In detail, B: B ~> (C is the unique linear map 
sending Lq^^DI G B to /(^) G (D; confer Lemma[5]for the uniqueness statement. 
In the sequel these identifications will be implicit. 

For future reference, we mention also that the usual Taylor polynomial allows 
one to provide a natural interpolator for the initial data, namely 

S*(/i, . . . , /™) = fi{x) + t f2{x) + ■■■+ (£^ /™(x), (3) 

which we will not need here because compute the kernel projector directly from 
its first-order factors. 

In this paper, we will study the Cauchy problem ^ for a completely reducible 
operator T{D), meaning one whose characteristic polynomial T(A) G ^[A] = 
(D[Ai, . . . , A„] splits into linear factors. Hence assume T = 7^™i...j'™'= with 
first-order operators Ti , . . . , T^ G (D [D] . By a well-known consequence of the 
Ehrenpreis-Palamodov theorem, the general solution of Tu = is the sum of 
the general solutions of the factor equations T™^u = 0, . . . , Tjf''''u — 0; see the 
Corollary on [51 p. 187]. Hence it remains to consider differential operators that 
are powers of first-order ones (we may assume all nonconstant coefficients are 
nonzero since otherwise we reduce n after renaming variables) . 



Of course the problem may still be ill-posed; we will not treat this issue here. 
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Lemma 13. Let T = a + a^Dt + aiDi + • ■ • + anDn e ^[D] be a first-order 
operator with all ai ^ 0. Order the variables such that all cumulative sums oq + 
Oi + • ■ • + fli-i are nonzero. Then the general solution ofT"^u — is given by 

m 
U[t, Xi, . . . , Xn) ^ / ^ Ci[Xi, . . . , Xn) (i — l)\ ' \ •' 

Xi =t + xi-\ h a;j_i - (ao + ai H h aj_i) xjai, (5) 

where (/i, . . . , /„) £ -7^,"Li are arbitrary functions of the indicated arguments. 

Proof. This can be found in some textbooks on differential equations [1, p. 139]. 
Setting up the change of variables given by (O and i = t, the differential operator 
is r = {a + ao DiT in th^ new coordinates (the ordering ensures invertibility) . 
Clearly, its fundamental solutions are u{f) = q ?~^e~''*/"° ii = 1,...,to), and 
the "integration constants" Ci are arbitrary functions of xi, . . . , a;„. D 

In principle, one could now combine the general solutions (0]) for each factor, 
substitute them into the initial conditions of ^ and then solve for the Ci in 
terms of the prescribed boundary data (/i, . . . , /„). With this choice of c^, the 
general solution will become the state operator for the Cauchy problem. However, 
this is a very laborious procedure, and therefore we prefer to use another route. 
Since we assume a completely reducible operator, we can employ the product 
representation of Proposition [HI In that case, it remains to consider the case of 
a single first-order factor. 

Lemma 14. Let T = a + a^Dt + aiDi + • ■ • + anDn G '^[D] be a first-order 
operator with all at ^ 0. Then the Cauchy problem Tu — 0, u(0,xi, . . . ,x„) = 
/(xi, . . . ,a;„) has the state operator H{f) — 6""*'"° Z* . Moreover, the differen- 
tial operator T has the right inverse 

T* = a-^ e°*/"« Z* At e""*/"" Z* . 

Here Z £ C"^]^ is the transformation ([5]) with t — t, and Z is its inverse. 

Proof. This follows immediately from Lemma [T3l The right inverse is computed 
using Lemma 3 of 11, after transforming the LPDE to a LODE. D 

By Propositio rilli we can determine the kernel projector for the Cauchy 
problem of Lemma [TJ] as P = H o trc, where trc('u) = m(0, xi, . . . , x„) in this 
simple case. Having the kernel projector and the right inverse T* in Lemma [T^ 
the signal operator is computed by G = (1 — P)T^ as usual. Now we can tackle 
the general Cauchy problem ([2]) by a simple special case of Proposition [HI 

Proposition 15. Let Ti,T2 € C[-D] be two first-order operators with nonzero 
coefficients for Dt. If Lq^^ is the evaluation defined after Theorem. \1'A then we 
have 

(Ti, [io.e I e G IR]) (7^2, [io,e I e e K]) = (riT2, [Lo,6^o.«A \^^n 

for the product of the Cauchy problems. 



12 



Proof. By the definition of tlie product of boundary problems, we fiave to show 
that [LQ,^,Lo^^Dt] — [Lq^Ci -^0,5^2]. Since each of these is defined as a biortho- 
gonal, it suffices to prove that the system m(0,^) = ut{0,£,) = has the same 
solutions as the system m(0,^) = (r2M)(0,^) = 0. But the latter is given by 

where the first term and the x^-derivatives vanish since u(0, ^) = 0. Using gq ^ 0, 
this implies that the two systems are indeed equivalent. D 

This settles the completely reducible case: Using Proposition[T51 we can break 
down the general Cauchy problem (U) into first-order factors with single ini- 
tial conditions. For each of these we compute the state and signal operator via 
Lemma [Ml hence the state and signal operator of (O by Proposition [Q] 



5 Conclusion 

As explained in the Introduction, we see the framework developed in this paper 
as the first stage of a more ambitious endeavor aimed at boundary problems for 
general constant-coefficient (and other) LPDEs. Following the enumeration of 
the Introduction, the next steps are as follows: 

1. Stage (1) was presented in this paper, but the detailed implementation for 
some of the methods explained here is still ongoing. The crucial feature of 
this stage is that it allows us to stay within the (rather narrow) confines of 
the PIDOS algebra. In particular, no Fourier transformations are needed in 
this case, so the analytic setting is entirely sufficient. 

2. As we enter Stage (2), it appears to be necessary to employ stronger tools. 
The most popular choice is certainly the framework of Fourier transforms 
(and the related Laplace transforms). While this can be algebraized in a 
manner completely analogous to the PIDOS algebra, the issue of choosing 
the right function space becomes more pressing: Clearly one has to leave 
the holomorphic setting for more analysis-flavoured spaces like the Schwartz 
class or functions with compact support. (As of now we stop short of using 
distributions since that would necessitate a more radical departure, forcing 
us to give up rings in favor of modules.) 

3. For the treatment of genuine boundary problems in Stage (3) our plan is to 
use a powerful generalization of the Fourier transformation — the Ehrenpreis- 
Palamodov integral representation [2], also applicable to systems of LPDEs. 

Much of this is still far away. But the general algebraic framework for bound- 
ary problems from Sections [5] and |3] is applicable, so the main work ahead of us 
is to identify reasonable classes of LPDEs and boundary problems that admit a 
symbolic treatment of one sort or another. 
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